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Through variational Monte Carlo calculation on Gutzwiller projected wave functions, we study 
the quasiparticle(qp) weight for adding and removing an electron from a high temperature super- 
conductor. We find the qp weight is particle-hole symmetric at sufficiently low energy. We propose 
to use the tunneling asymmetry as a tool to study the mechanism of electron incoherence in high 
temperature superconductors. 

PACS numbers: 



The scanning tunneling niicroscopy(STM) plays an im- 
portant role in the study of the high temperature super- 
conductors since it can provide local information on the 
single particle properties with ultrahigh energy resolu- 
tion. A striking feature in the STM spectrum of the 
high temperature superconductors is their remarkable 
particle-hole asymmetry. The hole side of the spectrum 
always dominate the particle side of the spectrum in hole 
doped cuprates 1]. 

This asymmetry is not at all surprising if we take the 
high temperature superconductors as a doped Mott in- 
sulators described by the t — J model. In such a doped 
Mott insulator, an added electron has a reduced prob- 
ability to contribute to the electron spectral weight in 
the low energy subspace of no doubly occupancy. More 
specifically, if the hole density in the system is x, then the 
total spectral weight in the particle side of the spectrum 
is reduced to x, while the total spectral weight in the 
hole side of the spectrum is not affected by the no dou- 
ble occupancy constraint. Thus the total spectral weight 
is particle-hole asymmetric for small x. However, such 
an asymmetry on the total spectral weight tell us noth- 
ing about the distribution of the spectral weight at low 
energy. To address the problem of tunneling asymmetry 
in the near vicinity of the chemical potential, we need 
more detailed information on the low energy excitation 
of the system. 

Rantner and Wen addressed this problem in the slave- 
Boson mean field theory of the t ~ J model 2]. In their 
theory, the tunneling asymmetry comes from the inco- 
herent part of the electron spectrum. In the slave-Boson 
mean field theory, an electron is split into two parts, 
namely the Fermionic spinon part that carry spin and 
the Bosonic holon part that carry charge. The supercon- 
ducting state is described by the Bose condensation of 
holon in the background of BCS pairing of spinon. In 
the presence of the holon condensate, the electron spec- 
trum acquire a nonzero quasiparticle{qp) weight. In the 
mean field theory, the spectral weight in the particle side 



of the spectrum is totally coherent since the holon re- 
moved during the particle injection process must origi- 
nate from the holon condensate. However, the hole side 
of the spectrum involves both coherent and incoherent 
part since the holon injected during the particle remov- 
ing process can stay either in the condensate or out of it. 
In the slave-Boson mean field theory, the qp weight dif- 
fers from the BCS result by a constant renormalization 
factor X, namely xu'^ for adding an electron and xv'^ for 
removing an electron. Thus if one neglect the asymmetry 
due to the band structure effect near the chemical poten- 
tial, then the qp weight is particle-hole symmetric and 
tunneling asymmetry must originate from the incoherent 
part of the electron spectrum. 

Recently, Anderson and Ong addressed the same prob- 
lem with a variational approach 3J. They constructed the 
variational wave function for the ground state and the ex- 
cited state of the t — J model following the original RVB 
idea. In their treatment, the particle-hole asymmetry is 
taken into account explicitly in the variational wave func- 
tion by the introduction of a fugacity factor Z(= ]^)- 
This fugacity factor is expected to play the role of the 
Gutzwiller projection into the subspace of no double oc- 
cupancy. In this theory, the electron spectrum is dom- 
inated by the qp contribution and the tunneling asym- 
metry originates from the coherent rather than the inco- 
herent part of the electron spectrum. Especially, when 
the excitation energy considered is much larger than the 
pairing gap (so that the particle-hole mixing due to the 
superconducting pairing is negligible), the qp weight for 
adding an electron is a factor Z smaller than that for 
removing an electron. In this theory, the qp weight for 
adding an electron scales linearly with the hole density x 
near half filling , while the qp weight for removing an elec- 
tron scales linearly with 1~ x. The problem of tunneling 
asymmetry and the qp weight is also addressed variation- 
ally in some other recent works |j,|3, la tD]- For example, 
it is proved by Yunoki that the electron spectral weight 
in the particle side is exhausted by the qp contribution 



for a superconductor described by Gutzwiller projected 
BCS wave functions J5|- However, a clear understanding 
of the hole-like qp excitation of the Gutzwiller projected 
state is still absent. 

In this paper, we show through the variational Monte 
Carlo calculation on Gutzwiller projected wave functions 
that the qp weight in the t-J model is particle-hole sym- 
metric at sufficiently low energy. Especially, we find the 
qp weight for adding and removing an electron in a pro- 
jected Fermi sea state both converge to the value de- 
termined from the jump of the momentum distribution 
function on the Fermi surface. We find the qp weight 
for adding and removing an electron both vanish at zero 
doping, as predicted by the slave-Boson mean field the- 
ory. However, the qp weight calculated from projected 
wave function show more complex momentum and dop- 
ing dependence. We find the qp weight vanishes like x"^ 
near half filling. The exponent x"'' increases monoton- 
ically from i to 1 when one move from deep inside the 
Fermi surface to far outside it in two spatial dimension. 
Thus the qp weight vanishes more slowly near zero doping 
than that predicted by the slave-Boson mean field the- 
ory. According to our result, the tunneling asymmetry 
near the Fermi level should be attributed to the incoher- 
ent part of the electron spectrum. We propose to use the 
tunneling asymmetry to study the mechanism of electron 
incoherence in the high temperature superconductors. 

In the Landau theory of Fermi liquid, the quasiparticle 
plays a dual role. On the one hand, the qp can be thought 
of as a particle-like elementary excitation on the ground 
state of a A^ particle system. On the other hand, the qp 
can also be thought of as a constituent of the grounds 
state of the N + 1 particle system, provided that the qp 
is located on the Fermi surface. The qp weight for adding 
an electron into the system on the Fermi surface is thus 
equal to the square of the matrix element of the electron 
creation operator between the ground state of N particle 
system and the ground state of TV -I- 1 particle system. 
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while the qp weight for removing an electron from the 
system on the Fermi surface is equal to the square of the 
matrix element of electron annihilation operator between 
the ground state of N particle system and the ground 
state of A — 1 particle system. 

In the thermodynamic limit, we have 
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thus the qp weight should be particle-hole symmetric. 
This simple argument do not apply to the superconduct- 
ing state. However, we do not expect the superconduct- 
ing pairing to change the conclusion. The quasiparticle 



in the superconducting state is a mixture of particle and 
hole components. Thus the superconducting pairing is 
expected to enhance rather than reduce the particle-hole 
symmetry. 

Now we calculate the qp weight variationally with the 
Gutzwiller projected BCS-type wave functions. This 
kind of variational description is widely used in the study 
of the t-J model and is believed to capture the low en- 
ergy physics of the cuprates[8l|9|,|l3,lll|, especially when 
the doping is near optimal. The variational ground state, 
namely the Gutzwiller projected BCS state with N par- 
ticle is given by (in the following we follow the notations 
ofi) 
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(1) 



where Pn is projection onto the subspace of N parti- 
cles, Pq is the projection onto the subspace of no double 
occupancy. |BCS) denotes the unprojected BCS mean 
field ground state. The elementary excitation above the 
ground state can be similarly constructed by Gutzwiller 
projection of BCS mean field excited state. For example, 
an particle-like elementary excitation is constructed as 
follows 
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where •f^ ^ is the creation operator for Bogliubov quasi- 
particle on the BCS mean field ground state which is 
related to the original electron operator c^ ^ through the 
Bogliubov transformation 
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Similarly, a hole like elementary excitation with the same 
momentum and spin is constructed as follows 
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The qp weight for adding and removing an electron are 
given by the overlap between the corresponding bare elec- 
tronic states and the qp states 
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As pointed out by Yunoki, using the fact that 
PgC^^Pg = ^GCkCT01 ' the qp weight for adding an 
electron can be related to the momentum distribution 
function Uk in the following way. 
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where ^ denotes the mean value of the particle num- 
ber with opposite spin. Thus to calculate Z^ , one only 
need to evaluate an ground state expectation value. The 
calculation of Z^ is more complex. However, using 
PqcI ^Pq — PgcI „ and the Bogliubov transformation 
Eq.(4), we are able to show that 



Z^=a 
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where the constant a is given by 
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and plays the role of the fugacity factor. Ok is given by 
following overlap integral 



Ok = 



(VI/N|VI/N) 



in which 
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Further simplification is possible when there is no su- 
perconducting pairing. In this case, one find 



where Arikp denotes the jump of Uk on the Fermi surface. 
Thus both Z^ and Z^ converge to Arikp on the Fermi 
surface in the thermodynamic limit, consistent with a 
Fermi liquid interpretation of the quasiparticle excita- 
tion. 

We now present the result of VMC calculation. Fig.(l) 
shows the momentum distribution function and the qp 
weight for the Gutzwiller projected Fermi sea on a 18 x 18 
lattice. Here the mean field Fermi sea is generated by 
a nearest neighboring hopping term on a square lat- 
tice with periodic-periodic boundary condition. The hole 
number is kept at 42(a; ~ 0.13) so that the closed-shell 
condition is satisfied. As shown in Fig.l, the total qp 
weight is a continues function of momentum across the 
Fermi surface. Thus the qp weight must be particle-hole 
symmetric on the Fermi surface in the thermodynamic 
limit, as we have argued before(note that the qp weight 
for adding (removing) an electron vanishes inside (outside) 
the Fermi surface). Note also that at this doping level 
the qp weight for removing an electron is only slightly 
higher than that for adding an electron in the whole 
Brillouin zone. Thus the quasiparticle contribution to 
the STM spectrum should be particle-hole symmetric at 
energy scale much smaller the band width. 
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for k inside the Fermi surface. Noticing the fact that 



N+i\ ('ivj/^ 1^) is nothing but the variational ground 
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state of the N -I- 1(A^ — 1) system for k located on the 
Fermi surface , we have 
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where kp denotes momentum located on the Fermi sur- 
face. Thus, if the qp weight is a continuous function of 
particle number, it should be particle- hole symmetric in 
the thermodynamic limit. Furthermore, using Yunoki's 
relation and the fact that Z^ vanish for k inside the Fermi 
surface, we have 



FIG. 1: VMC result for Gutzwiller projected Fermi sea on a 
18 X 18 lattice with 42 holes(2; ~ 0.13). The hole density is 
choosen that the closed shell condition is satisfied for system 
with periodic-periodic boundary condition. The mean field 
state is generated by nearest-neighbouring hopping term on 
the lattice, (a) momentum distribution function, (b)total 
qp weight. (c)qp weight in the (0,0) — (7r,7r) direction. (d)qp 
weight in the (0, 0) — (0, vr) direction. 
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for k inside the Fermi surface and 
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Another implication of the above result is that most 
spectral weight in the hole side of the spectrum is inco- 
herent. To show this more clearly, we plot in Fig. (2) the 
doping dependence of the -^1^(00) ^^^ ^k^fiTTr)- F^o™ 
the figure we see that both ■^k=(o o) ^"^^ '^\=(t: tt) ^^'^" 



ish near half filling, in agreement with the slave-Boson 
mean field theory prediction. As compared to the slave- 



Boson mean field result ( Z^ 



- ^,,2 7- 



'^k = 2;wj), the 



qp weight calculated from projected wave function show 
more complex momentum dependence. The origin of the 
momentum dependence can be traced back to the non- 
monotonic behavior of n{k). As first noted by in |l2l| . 
the non-monotonic behavior of n(k) comes from the cor- 
related hoping nature of the electron in the t-J model. 
Another difference between the slave-Boson mean field 
theory and the projected wave function is the doping de- 
pendence of the qp weight (although they both vanish at 
half filling) . In the slave-Boson mean field theory, the qp 
weight scales linearly with x at all momentums. How- 
ever, we find the qp weight calculated from the projected 
wave function show x"^ behavior at low doping, where 
Qfk is a momentum dependent exponent. For the two 
dimensional projected Fermi sea, ak increase monotoni- 
cally from i deep inside the Fermi surface to 1 far out- 
side the Fermi surface. Thus the qp is more robust in the 
variational description than in the slave-Boson mean field 
theory. This difference may be caused by spin-charge re- 
combination induced by the Gutzwiller proiectionjl3j. 
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FIG. 2: Doping dependence of the qp weight as calculated 
from the Gutzwiller projected Fermi sea on a 18 x 18 lattice 
as a function of hole density. Periodic-antiperiodic boundary 
condition is used in the calculation to maximize the number of 
hole density that satisfy the closed-shell condition, (a)doping 
dependence of .Zk at three momentums along (0,0) to (7r,7r), 
namely k = (0, 0),kFfflnd(7r,7r)(note that since the boundary 
condition is periodic-antiperiodic, the momentums are not ex- 
actly in the nodal direction). (b) momentum dependence of the 
exponent Ok in the nodal direction. 



Next we present the result for the Gutzwiller projected 
d-wave BCS state. The mean field state is generated by 
the following Hamiltonian 
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in which d^ is the form factor for d-wave pairing and 
the sum is limited to nearest-neighboring sites. We take 
A = 0.1 and the hole number is kept at 42. \i is de- 
termined by the mean field equation for density. Fig. (3) 
shows the momentum distribution function and the qp 
weight on a 18 X 18 lattice with periodic-periodic bound- 
ary condition. The result is basically the same as that 
of the projected Fermi sea state apart from the particle- 
hole mixture in the vicinity of the Fermi surface. One 
thing to note is that the total qp weight is a monotonic 
function of momentum in the nodal direction and do not 
exhibit pocket structure of the kind found for the Zy in 
[6j. To clarify the situation, we plot the momentum de- 
pendence of Z^ and Z^ separately in Fig. (4). While a 
well defined pocket structure is observed in Zy , no corre- 
sponding structure exist in Z^ . The origin of the pocket 
structure in Z^ can be traced back to the non-monotonic 
behavior of n(k) around the nodal point as a result of the 
correlated hopping of electron in the t-J model and the 
d-wave pairing. However, Z^ do not show such pocket 
structure since it is a monotonically decreasing function 
of momentum in the nodal direction. At the same time, 
since the decrease of Z^ in the (0,7r) to (7r,0) direction 
due to pairing is compensated by the increase of Z^^ , the 
total qp weight do not show pocket structure. 
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FIG. 3: VMC result for Gutzwiller projected d-wave BCS 
state on a 18 X 18 lattice with 42 holes(a; ~ 0.13) and y = 0.1. 
(a) momentum distribution function, (b)total qp weight. (c)qp 
weight in the (0,0) — (ti", tt) direction. (d)qp weight in the 
(0, 0) - (0, tt) direction. 

Above we have show that qp weight is particle-hole 
symmetric at sufficiently low energy. Thus the tunnel- 
ing asymmetry at low energy should be attributed to the 
incoherent part of the electron spectrum. Since the elec- 
tronic spectrum in the particle side is totally coherent, by 



tunneling asymmetry at low energy. 
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FIG. 4: Momentum dependence of Z^ and Zy. in the 
Gutzwiller projected d-wave BCS state. 



subtracting the STM spectrum on the hole side from that 
on the particle side we can extract the incoherent part 
of the electron spectrum. For strongly correlated system 
like cuprates, the information on the incoherent spectral 
weight is of great value . Through studying this spec- 
trum, one can figure out the mechanism by which the 
bare electron decay into the many particle excitations 
and also the nature of the many particle excitations it- 
self. In the context of the cuprates, two hotly discussed 
mechanisms to generate electron incoherence are scatter- 
ing with some bosonic collective mode(like the neutron 
resonance mode) |l4| and electron fractionalization as we 
have discussed in the slave-Boson mean field theory |2|. 
It is interesting to see which mechanism dominate the 
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